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I. Tasks Performed 

The major project carried out under this grant was the 
measurement of electric current on two Nike Tomahawk sounding 
rockets launched by the Swedish Space Corporation from ESRA'NGE, 
Kiruna, Sweden on 7 February 1972 and 9 February 1973. The 
instruments used to make the measurement were split Langmuir 
probes. The work done may be divided into three phases: 
experiment design, construction and flight, and data analysis. 

The experiment design work included continuing efforts to under- 
stand the performance and behavior of the split Langmuir probe 
by analyzing data obtained from previous flights. The scientific 
results of this work, including previous flight results only partially 
supported by this grant, are summarized in section II of this report. 
The details of the scientific results of the Swedish flights alone 
are given in section III. The details of the instrument design 
which was developed for this project are given in section IV. 

Most of the work done under this grant was carried out by 
a graduate student, Edgar Bering, who used the results for his 
Ph. D. thesis entitled "Auroral Zone Electric Current Measurements 
on Sounding Rockets" (University of California, Berkeley, 1974). 

II. Summary of Scientific Results 

An instrument, the split Langmuir probe, has been developed 
to make in situ measurements of ionospheric current density and 
plasma bulk flow. The split Langmuir probe consists of two 
conducting elements that are separated by a thin insulator that 
shield each other over a 2rr solid angle, and that are simultan- 
eously swept from negative to positive with respect to the plasma. 

By measuring the current to each plate and the difference current 
between plates, information can be obtained on the plasma's 
current density, bulk flow, electron temperature, and density. 
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Sources of error limiting complete measurement of all of these 
parameters include plate area differences, plate work function 
differences, input electronics differences, and probe wake effects. 
The instrument was successfully flown twice on sounding rockets 
into auroral events near local midnight on Swedish Nike Tomahawk 
7/1 from ESRANGE, Kiruna, Sweden on 7 February 1972, and on 
Swedish Nike Tomahawk 7/2 from ESRANGE on 9 February 1973. 

In addition, data obtained from the split Langmuir probes flown 
on Javelin 8. 56 VE from Fort Churchill, Manitoba, Canada on 
3 April 1970 was analyzed. On the Javelin flight, the rocket 

passed repeatedly in and out of a bright folded arc structure. In 

the southern edge of the arc, the parallel detector measured a 

- 4 2 

downward thermal ion current of 1.7x10 amperes/m parallel 
to B in the earth-fixed frame. The latitudinal extent of this 
enhanced parallel flow was a few hundred meters. The upward 
return current was carried by precipitating electron within and 
to the north of the visual form. Current balance implies that 
parallel electron mobility was reduced in the region of the ion 
beam by an order of magnitude relative to expectation, apparently 
by anomalous resistivity.- On the first Swedish flight, the electric 
field deduced from the ion flow velocity measurements and the 
electric field measured by the double probe experiment agree 
at altitudes above 140 kilometers to an accuracy well within 
the uncertainties of the two measurements. Below 140 kilometers, 
the two measurements diverge. This divergence gave a measure- 
ment of perpendicular current density, conductivity, and neutral 
winds. The height integrated current was 0.11 amperes/m flowing 
at an azimuth of 276°. This current was an order of magnitude 
too small to have been the main electro jet. The neutral winds 
were strong and to the south. On the second Swedish flight, 
the current measured on the upleg was 0. 68 amperes/m flowing 
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at an azimuth of 55°. On the downleg, an equally strong south- 
westward current was indicated, implying the presence of oppo- 
sitely directed electrojets. The neutral winds were strong and 
highly structured. The instrument is also ideal for studying the 
near wake of a small probe. Significant wake perturbations in 
plasma density and temperature were found due to the probe body 
itself, even though the probes were on the order of or smaller than 
the Debye length. The wake of small probes show apparent mag- 
netic field aligned structure, even though the probes were much 
smaller than the ion gyroradius. 

III. Results From Swedish Nike Tomahawks 7/1 and 7/2 

A. Introduction 

1. General Comments 

The most recent two flights of the split Langmuir probe 
were on two identical Nike -Tomahawk payloads launched from 
ESRANGE, Kiruna, Sweden near local midnight on 7 February, 
1972 and 9 February, 1973. These rockets were the seventh 
Nike Tomahawk payload configurations to be launched by the 
Swedish space research program and are customarily referred 
to as SNT 7/1 and 7/2. 

A sketch of the payload is shown in Figure 1. Probe 
geometry and sweep characteristics are given in Tables I and II. 
Details of the probe and the instrument electronics are dis- 
cussed in section IV. The other experiments on the payload 
relevant to auroral physics included a double probe electric 
field experiment provided by the Royal Institute of Technology, 
Stockholm; Langmuir probes, photometers, and a Faraday 
rotation experiment from the Uppsala Ionospheric Observatory; 
and an array of low energy electron and proton detectors from 
the Kiruna Geophysical Observatory. 
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2. SNT 7/1 Flight Conditions and Payload Performance 

Payload 7/1 was launched into the recovery phase of a 
400 v auroral substorm. Figure 2 shows the ground magne- 
tometer record from the period around the flight. A launcher 
malfunction sent the rocket on a course that was nearly due 
west. It is somewhat difficult to determine where the aurora 
was in relation to the rocket because darkroom mishaps 
destroyed the all-sty camera records from the sites that were 
clear, and there was broken cloud at the sites from which 
the various investigators were observing. Stoffregen (private 
communication) reported the auroral forms as corona. As 
viewed through broken cloud from the launch site, the aurora 
appeared to be brighter to the north of the trajectory, and to 
be brighter in the west than the east. It is hoped that analysis 
of the particle and photometer data will provide more infor- 
mation on the relative location of the rocket and the aurora 
than is now available. 

There were several problems which limited, but did 
not seriously impair, experiment performance on the SNT 
7/1 flight. As the result of an error in interpreting payload 
potential results from previous flights, the sweep stepped 
over plasma potential instead of sweeping through it. There- 
fore, it was impossible to measure electron temperature and 
plasma density or to attempt to make a direct measurement of 
current density. Another sweep potential related problem 
arose because the electron current being drawn from the plasma 
by the various experiments was rather large. Asa result, 
the payload potential was very negative for about 40 seconds 
early in the flight. This effect interacted with the reference 
probe low pass filter described in section IV in such a way 
that the sweep was always positive from 80 to 120 seconds, 
and no ion flow data was obtained for this time period. 
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Payload mechanical performance also created a problem 
for the split probe experiment. The spin rate acheived was 
higher than nominal. The parallel probe was motor 
mounted, and was supposed to rotate 180° every 4 sweeps in 
order to measure any DC offsets. Because of weight con- 
straints, the motor used generated only marginally enough 
torque to move the probe through a potential minimum at 90° 

rotation caused by the probe's flexing. At the higher spin 
rate, the motor torque proved insufficient, and the parallel 
probe stuck 90° out of alignment. 

3. SNT 7/2 Flight Conditions and Payload Performance 

Payload 7/2 was launched into a band of at least three 

parallel quiet arcs (Figure 3). These arcs occured during 

a relatively quiescent interval between two substorms on a 

moderately active night. The magnetogram from ESRANGE for 

that night appears in Figure 4. K was 4+ for the three hours 

P 

before midnight and 4- for the three hours after midnight. 

Figure 4 indicates that there was about a 200y magnetic 
disturbance at the time of launch. The rocket trajectory 
went almost north, and the projection for the trajectory down 
to 205 km along B traverses the southernmost half of the more 
intense visual forms appearing in Figure 3. 

Malfunctions severely limited the quality and quantity 
of the data returned by the split Langmuir probe on SNT 7/2. 

The major problem was the appearance of a substantial 20kHz 
noise band on the PCM telemetry channel. The noise appeared 
after discrimination, and was therefore added to the signal 
before it entered the telemetry electronics. What appears 
to be the same noise band was detected by a swept frequency 
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analyzer measuring the AC electric field. From the character 
of the signal, it appears that it is being radiated by the rocket. 

It is believed that the cause of this noise was a time variable 
resistive break in the signal ground connection between the 
split Langmuir probe and telemetry. The effect of the noise 
on the split Langmuir probe data was to make the signal quite 
difficult to read after t = 120 seconds, and useless after 
t = 130 seconds except for a six second burst at t = 350 seconds.' 

Another significant problem was that the total plate to 
plate offset voltage was greater than 0. 1 kT^/e. This caused 
the DC offset in the difference current at plasma potential to 
be outside the dynamic range of the telemetry system, so that 
it was again impossible to make a direct measurement of 
current density. 

B. Single Plate Langmuir Probe Measurements 

No single plate data analyses were performed on the SNT 7 
data. For the first flight, this was impossible since the sweep 
did not pass through plasma potential. On the second flight, 
the effect of sweeping slowly with respect to the rocket spin 
was to introduce a spin modulation on the single plate data 
that made analysis exremely difficult. The plasma density 
was also measured on this payload by several experiments 
provided by the Uppsala Ionospheric Observatory (UJO). 

Because of the difficulties with the split probe, the UJO data 
have been used for all of the analyses described in this chapter. 
The UJO results are plotted as a function of altitude in Figures 
5 and 12., This data was provided with remarkable speed 
through the courtesy and efficiency of Bengt Holback. 
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C. Difference Current Measurements at Large Negative 
Potentials -- Ion Flow Measurements: Results From 
SNT 7/1 

1. Methods of Data Analysis 

The difference current data from SNT 7/1 were analyzed 
using the methods described in Bering (1974) with one addition and 
one exception. The addition to the analysis procedure resulted from 
the fact that the quasi-cylindrical geometry employed behaved very 
much like a cylinder in its focusing properties. Therefore, the 
focusing factor was modeled using the results of Bettinger and 
Walker (1965) with a single point normalization as described in 
Chapters III and V of Bering (1974). The exception arose when 
it became apparent that the wake correction described in section 
III. A. 2. f of Bering (1974) increased the discrepancy between the 
split and double probe data rather than decreasing it. This cor- 
rection is only used to illustrate the possible range of systematic 
errors. 


2. Parallel Flow 

The inability of the parallel probe motor to turn the parallel 
probe means that only the two components of ion flow velocity 
normal to the spin axis of the rocket have actually been 
measured. In order to transform these data into the earth- 
fixed geomagnetic reference frame it is necessary to determine 
the flow velocity parallel to the rocket. Since the rocket 
precession period was 22 seconds, it is possible to use the 
precession analysis method described by Bering. The 
method involves deducing the perpendicular components of the 
apparent electric field in earth -fixed geomagnetic coordinates 
for various assumptions of v^, the component of plasma velocity 
parallel to B. The correct choice of v g is the one which mini- 
mizes the apparent precession periodicity in the deduced 
electric field. 
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The electric field components deduced from the SNT 7/1 
ion flow data for various assumptions of Vg appear in Figure 6. 
Positive Vg is upward along B. Detailed study of the results 
presented in Figure 6 indicates that a choice of v D that is 
time varying and substantially different from zero is required 
to minimize apparent precession periodicity in the results. 

The results of this study appear in the bottom panel of Figure 7 

labeled "best fit" and in the bottom panel of Figure 8. 

The values range from +1 km/sec before t = 150 seconds to 
-0. 6 km/sec around t = 320 seconds, and are somewhat 
antisymmetric about apogee. 

These apparent parallel flow velocities present a sub- 
stantial problem in interpretation. The maximum apparent 
parallel ion velocity, 1 km/sec, is an order of magnitude greater 
than is seen by the Chatanika radar (Baron, private communi- 
cation; Mozer and Cullers, private communication). The 
anti-symmetry about apogee is also very suspicious. 

There are four possible interpretations of the v^ results. 
First, they could represent a real geophysical situation. This 
is extremely unlikely in view of the Chatanika results and the 
apogee anti- symmetry . Second, if the data before t = 160 
seconds is dismissed as the result of outgassing and other 
post -deployment settling down effects (Haerendel et al . , 1974), 
the data is not outside the realm of what has been seen before, 
and could represent a real geophysical situation. The mag- 
nitudes shown are still large enough for this interpretation 
to be regarded as substantially unlikely. Third, this result 
could be the result of a subtle malfunction internal to the 
instrument. The excellent agreement between the split and 
double probe data reported below argues strongly against this 
interpretation. Fourth, it could be caused by a rocket- 
plasma interaction of some sort which perturbs the component 
of ion velocitv Darallel to B. 
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As can be seen in Figure 7, the best fit results for 
are approximately consistent with the ion velocity parallel 
to B being zero in the rocket frame prior to t = 300 seconds. 
Similar results have been seen by another experimentor using 

a different technique on other flights (M. C. Kelley and J. -J. 
Berthelier, private communication). This other technique 

involves deducing the ion flow velocity from the phase of the 
wake observed with an ordinary probe at three positions around 
the precession cone of the rocket. Analyses of sounding rocket 
probe and particle data using this technique have given a similar 
result- -namely , that the parallel ion velocity is zero in the 
rocket frame, or nearly so. Results from this technique do 
not always agree with double probe results from the same rockets, 
and must, therefore, be treated with extreme caution. 

Nonetheless, it does appear that the fourth interpretation of 
the v^ results--that they are an artifact of the rocket-plasma 
interaction--is the most likely one. 

This phenomenon has not been previously reported in 
the literature. In view of its extreme interest to many sounding 
rocket investigators, it deserves more investigation. Several 
other properties can be deduced from the available split 
Langmuir probe data. First, it is not a universal phenomenon. 

The split Langmuir probe data gives no evidence of the effect 
on the Black Brant (Bering et al . , 1973) or Javelin rockets (Bering, 
1974). Second, it appears to be a DC phenomenon, and to be a 
slowly varying function of rocket velocity. It certainly does not 
prevent abrupt changes in v B from being detected, and probably 
does not affect the validity of results based on detection of such 
changes in parallel difference current. Finally, it does not appear 
to be related to or to invalidate measurements of parallel electric 
fields, since the effect has been observed both on. flights where 



11 . 


large E|| was observed and on flights where large E|| was not 
observed. Conversely, it has been absent in both circumstances. 
In conclusion, the SNT 7/1 parallel flow results appear to indicate 
the presence of a new and puzzling rocket-plasma interaction 
requiring further study before it can be fully understood. 

3. Perpendicular Flow 

The problem of interpreting the parallel results has 

implications as far as the perpendicular flow results are 

concerned. It is unclear which is the best procedure to follow 

in analyzing the data: use the best fit values for v , use 

v = v , or use v = 0. For purposes of comparison, 

B rocket, B B 

all three possibilities are plotted in Figure 7 along with the 
results of applying the wake correction to the best fit results. 

The major conclusion that can be drawn from Figure 7 is that 
it does not make an enormous amount of difference which 
analysis procediire is chosen. Therefore, the results using 
the best fit values for v,, will be used for the rest of the 

ID 

analyses described below. In order to illustrate the range 
of possible systematic errors in subsequent analyses using the 
perpendicular data, the analyses were performed at selected 
times on all the results shown in Figure 7.. 

The electric field was also measured on the SNT 7 
payloads by a double probe experiment provided by Ulf Fahleson 
of the Royal Institute of Technology, Stockholm, Sweden. Dr. 
Fahleson was kind enough to transmit the data from one pair 
of spheres via hardwire in the payload to the PCM encoder 
in the split Langmuir probe. In addition, he was kind enough 
to consent to analysis of this data at the University of California 
as part of this research, prior to his analyzing the data. 
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The double probe electric field data has been analyzed 
using the method of Kelley (1970) and compared with the results 
of the split Langmuir probe analysis (Figure 8). Above 140 
km in altitude, the two sets of results are in excellent agreement, 
being well within the errors in the measurements. This 
agreement is substantially better than any reported in previous 

chapters, and represents further confirmation that the two 
experiments are measuring the ambient electric field. The 

difference between the two results which develops below 140 
km on the downleg does not appear to be the result of errors 
in the measurement. Instead, it is due to the effect of ion- 
neutral collisions on ion mobility, and can be transformed into 
a measurement of perpendicular current density. 

4. Current Density 

The current density in an ionized medium can be deduced 
from the relationship 

j = ne ( v. - v ) (1) 

if the ion and electron velocities are known. In the ionosphere, 
electron mobility is not substantially affected by collisions 
above about 90 km. It is possible, therefore, to deduce the 
electron flow velocity above 90 km from the double probe 
electric field measurements using 

v = E X B / B 2 (2) 

e 

The SNT 7/1 ion flow and electric field data have been analyzed 
using Equations 14 and 15 to yield a measurement of perpendi- 
cular current density. The results of these calculations are 
plotted as a function of altitude in Figure 9. 
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The data indicates that the rocket passed through a pre- 
dominately westward electrojet, with a maximum current 
_ 6 2 

density of 5 x 10 amperes /meter oc curing at an altitude 
of 110 kilometers. The total height integrated current was 
0. 11 amperes /me ter flowing at an azimuth of 276°. The 
error bars shown in the figure represent the statistical 
uncertainty in the measurement. In order to estimate the possible 
impact of systematic errors, the data has been analyzed at 
selected points using the other possible values for the parallel 
flow and wake correction shown in Figure 7. The scatter in the 
results of the various analyses is not significantly greater than 
the calculated error bars. Comparison of the size of the error 
bars with the magnitude of the measured current leads to the 
conclusion that the split Langmuir probe experiment on SNT 7/1 
was accurate enough to make the current density measurement 
shown in Figure 9 valid and meaningful. 

The current density measurements are in qualitative agree- 
ment with previous work. For example. Park and Cloutier 
(1971) reported a magnetometer measurement of the electrojet 
in an aurora in which they deduce a current density of 1. 5 x 10 - ^ 
amperes /meter ^ in a westward electrojet. Also, they agree 
with theoretical electro jet models such as the one developed 
by Bostrdm (1964). 

Since substorms are highly variable phenomena and can 
vary in intensity over several orders of magnitude, it is far 
more useful to compare these results with the ground magneto- 
meter data from the same substorm. At t = 300 seconds, the 
magnetic disturbance at launch was Ax = -204, Ay = 116, Az = -60. 
The measured current is too small by an order of magnitude, 
and is oriented 3;5 too far to the north to have produced the 
measured ground disturbance, it appears, therefore, that 
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the rocket did not pass through the main electrojet. From the 
sign of z, and the observed location of the aurora, it appears 
that the rocket passed to the south of the region of maximum 
electro jet intensity. 

5. Conductivity 

Conductivity is formally defined by the relationship between 
electric field and current density. If both of these quantities 
are known, it is a simple task in algebra to calculate the 
conductivity. In a magneto-plasma, the problem is slightly 
complicated by the fact that conductivity is not a simple scalar, 
but rather a tensor. Perpendicular to B, there are two components 
of current density--parallel and perpendicular to the electric 
field- -and two different corresponding elements in the conductivity- 
tensor. These are known as the Pedersen and Hall conducti- 
vities respectively. These quantities have been calculated for 
the SNT 7/1 flight using the data shown in Figures 8 and 9. 

The results have been plotted as a function of altitude in Figure 

10 . 

The data indicate a maximum Hall conductivity of 1. 1 x 10 ^ 

mho/meter at 110 km and a maximum Pedersen conductivity of 
-4 

1.4x10 mho/meter at 117 km. The height integrated Hall 
conductivity was 3.0 mho and the height integrated Pedersen 
conductivity was 3.1 mho. The same comments about the size 
and relative importance of statistical and systematic errors 
that were made above in regard to current density apply here 
as well. It should be noted that the error bars are somewhat 
larger in percentage terms than in the previous section* a 
result to be expected since the measurement is an additional 
step away from the raw data. 

A model conductivity profile has also been plotted in each 
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panel of Figure 10 for purposes of comparison. Since the most 
recent related experiment to be reported is the one performed 
by Brekke et al . (1973), the model conductivity calculations 
were done wherever possible by the same methods in order 
to insure maximum comparibility of results. Specifically, the 
same approximation for the ion cyclotron frequency, the same 
ion-neutral collision rate coefficient, and the same neutral 
density model (Banks and Kockarts, 1973) were used. 

The comparison between the measured and computed 
conductivity- is, in general, quite good, and confirms the overall 
validity of the measurement. However, there are some dis- 
crepancies which are significant and have to be resolved. 

The measured Hall conductivity is somewhat smaller than predicted, 
the peak Pedersen conductivity is larger than predicted, and 
there are some alarming excursions into negative conductivity. 

The most likely explanation of these problems is that the effective 
conductivity in the earth-fixed frame is being changed by neutral 
winds. 


6. Neutral Wind 

The neutral wind velocity may be calculated from measure- 
ments of the ion flow velocity and electric field if the conductivity 
is known or assumed. In equation 2 of their paper, Brekke 
ejijal. give an expression for ion velocity involving neutral wind 
which may be solved for neutral wind as a function of ion flow 
and electric field. 


u 


v. 

i 


0 . 

<e+v XB> 
m 


(3) 


where u is the neutral wind velocity and Q. is the ion cyclotron 

frequency. Using the same values for £1. and v. that were used 

i in 
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above to compute the model conductivity profiles, the neutral 
wind velocity has been calculated from the SNT 7/1 data. This 
procedure is equivalent to assuming the conductivity to be equal 
to the calculated profiles shown in Figure 10. Figure 1 11 shows 
the results of these calculations plotted as a function of altitude. 

The neutral wind measurements indicate that the winds had 
somewhat high velocities and exhibited substantial altitude 
variation at the time of this flight. The north - south component 
rises from small values below 100 km to around 300 - 350 m/sec 
southward above 110 km. The west - east component varies 
substantially with altitude. It ranges from small values at 
high altitudes to 350 m/sec westward at 117 km and then 
reverses to 200 m/sec eastward at 105 km. At lower altitudes, 
the situation is obscured by rapidly increasing systematic 
errors, but the wind appears to become westward again. 
Generally similar comments to those made above can be made 
about the relative size and importance of statistical and sys- 
tematic errors. It should be noted that the relative size of the 
systematic errors appears to be somewhat larger in the neutral 
wind results than in the other results, and grows rapidly at the 
lower altitudes. The size of all errors grows rapidly above 
135 km. 

This is the first in situ point measurement of neutral winds 
ever made in the winter midnight auroral zone ionosphere 
making direct comparison with other data impossible. All of 
the other high altitude neutral wind measurement techniques 
either cannot be used in such circumstances, or do not make 
point measurements. The most commonly used technique, 
optical tracking of artificially introduced gas clouds (Smith, 

1968; Rees, 1971; S toff regen etal. , 1972; Stoffregen, 1972; 
Meriwether et al . , 1973), requires that the cloud be in sunlight, • 
a condition impossible to fulfill at winter midnight in the auroral 
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zone. None of the other three techniques--measurement of the 
Doppler shift in irglow emission lines (Hays and Roble, 1971), 
measurement of satellite drag (De. Vries, 1972) or incoherent 
scatter radar measurements of ion flow velocities in the E and 
F regions (Brekke et al . , 1973)--make point measurements. 

Despite the impossibility of a direct comparison, qualitative 
comparison of these results with the previous work cited above 
is still of interest. The range of magnitude and direction that 
have been reported at high latitude is enormous. The highest 
velocity measured, by De Vries, was 1000 m/sec. Roughly 
half of the results cited reported winds of 300 m/sec or more. 
There appeared to be a correlation of high velocities with magnetic 
activity, with proximity to the auroral zone, and with proximity 
to the midnight sector. Those authors who did show point 
measurement profiles from the auroral zone during aurora 
(Smith, Rees) reported profiles that were qualitatively similar 
to the results reported here in terms of variation with altitude. 

It seems, therefore, that the results shown in Figure 11 are 
in reasonable agreement with previous results and are a valid 
measurement of neutral wind velocity. 

D. Difference Current Measurements at Large Negative 
Potentials--Ion Flow Measurements: Results from 
SNT 7/2 

1. Methods of Data Analysis 

The methods used for analyzing the SNT 7/2 data were, 
in general, exactly the same as those used for SNT 7/1. There 
were some differences engendered by the noise problem discussed 
above. First, there was substantially more analog filtering used 
before passing the signal to the PCM decoder. Second, the 
criterion for synchronization used in preparing synch word 
checked tapes in the computer (Schutz, 1973) was relaxed 
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somewhat in order to facilitate correcting the timing on the 
final digital tapes. Finally, the sine wave fit procedure was 
monitored in detail in real time on a VISTA console, adjustments 
made in the control parameters as required, and spurious fits 
rejected from the output selected for final analysis. 

The analysis procedure was also altered somewhat by the 
fact that the parallel detector worked. The various DC offsets 
affecting the parallel detector, discussed in Chapter HI, were 
determined from the data by assuming that the ambient flow 
was constant or nearly so from one sweep to the next. The DC 
offsets were then solved for algebraicly at each probe input 
relay change and motor rotation. Since the assumption is 
somewhat questionable, this procedure is most accurate when 
the results are averaged over a large body of data centered 
on apogee. The noise problem made this impossible, and the 
accuracy of the resulting DC offset corrections is questionable. 

2. Parallel Flow 

Since the parallel probe data was not completely reliable, 
it was also necessary to determine the parallel ion velocity using 
the precession analysis method described above in Chapter IV, 
and in section C.l of this chapter. The electric field components 
deduced from the SNT 7/2 data for various assumptions of v 

B 

appear in Figure 13. The choice of v^ that is required to 
minimize the apparent precession periodicity in the results 
prior to t = 130 seconds is 1. 6 km/second. This result is 
plotted in the bottom panel of Figure 14 along with the velocity 
deduced from the parallel probe data, and with other possible 
assumptions as discussed above in regard to the 7/1 analysis. 
There is a substantial difference between the probe data and 
the results of the precession analysis. Because of the problems 
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with the DC offset corrections, the precession analysis result 
is believed to be more reliable. 

Comparison of the bottom panels of Figures 7 and 14 
indicate that the SNT 7/2 parallel flow data appears to be a 
more or less smooth continuation of the 7/1 results into the 
lower upleg region where no data was obtained on 7/1. This 
continuation is a remarkable coincidence, considering the 
difference in geophysical conditions. It is, in fact, no coin- 
cidence but rather a manifestation of the same rocket-plasma 
interaction effect discussed in section C. 2 of this chapter. 

3. Perpendicular Flow 

There are several problems one encounters in the perpen- 
dicular flow analysis that are difficult to resolve without the 
opportunity to compare the split and double probe results above 
140 km. They include choosing the best parallel flow analysis 
result, deciding whether or not to use the wake correction, 
checking the accuracy of the focusing model, and checking for 
any human errors in experiment installation. These possible 
problems, and the difficulty in resolving them, reduce the overall 
confidence in the results of the perpendicular flow measurement 
and subsequent analyses. Under the assumptions that no human 
blunders were made in preparing the payload, and that it is 
most appropriate to do so, the same analysis procedures and 
parameters have been emplqyed for SNT 7/2 as for 7/1. 

In addition to these methodological difficulties associated 
with the lack of data at apogee, the early upleg data from probe 
experiments may be subject to systematic errors of unknown 
magnitude duo to various "settling down" effects. These are 
hypothetical effects caused by outgassing and other adjustments 
to the space environment (Haerendel, et al . , 1974). The possi- 
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bility that these effects have perturbed the data also reduces con- 
fidence in the 7/2 perpendicular flow results. 

Despite these reservations, the quality of the noise-filtered 
7/2 data that does exist appears to be as good as the 7/1 data. 

No reason to doubt the validity of the 7/1 analysis choices as 
applied to the 7/2 data has arisen. All of the analyses have been 
carried through to completion, and the results merit serious 
consideration. 

The perpendicular flow data from the 7/2 flight appear in 
Figure 15, using the best fit precession analysis parallel velocity 
value as shown in the bottom panel. As above, the electric 
fields measured by Ulf Fahleson's double probe experiment 
and transmitted via the split Langmuir probe's PCM telemetiy 
system are also shown. Once again, there are differences between 
the two measurements which appear to be due to the effect 
of ion-neutral collisions on ion mobility and which can be trans- 
lated into a measurement of current density. 

4. Current Density 

The. current density has been determined from the SNT 7/2 
data using the same method as was used for the 7/1 data. The 
results are shown in Figure 16,. along with selected points from 
the other analyses shown in Figure 14. The same comments 
about the relative size and importance of statistical and sys- 
tematic errors can be made about these results as about the 
7/1 results, keeping in mind that the scatter in the results of 
the. various analyses might represent somewhat of an under- 
estimate of the systematic errors on this flight. 

On the upleg, the rocket passed through a north-east 

current, with a maximum northward current density of 1. 85 
-5 2 

x 10 amperes/meter at 112 km and a maximum eastward 
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current density of 4.4 x 10 amperes/meter at 105 km. 

The height integrated current up to 140 im was 0. 68 amperes/ 

meter, flowing at an azimuth of 55°. On the down leg, the single 

-5 

measurement made showed a current of 4. 2 x 10 amperes/ 
meter 2 flowing at an azimuth of 229°. These values are 
roughly comparable with previous measurements as cited in 
part C. 4 of this section. 

In comparison, the magnetic disturbance field measured 
at ESRANGE during the flight was Ax« -175y decreasing to 
-I25y, Ay fn 0, andAz« -70y. The measured magnetic distur- 
bance vector implied that the net ionospheric current was west- 
ward, with no north- south component. The ionospheric current 
actually measured by the rocket is not comparable to the current 
implied by the ground magnetometer data. This disagreement 
has some very interesting implications about auroral structure 
which will be discussed below. 

5. Conductivity 

The Hall and Pedersen conductivities have been calculated 
from the SNT 7/2 data using the same methods that were used 
for the 7/1 data. Figure 17 displays these conductivities, 
plotted as a function of altitude. Again, the same comments 
about errors apply. 

The maximum Hall conductivity on the upleg was 1.43 

_3 

x 10 mho/meter at 105 km. The maximum Pedersen conductivity 

_3 

was 1. 32 x 10 mho/meter at the same altitude. The height 
integrated Hall conductivity was 11. 7 mho and the height inte- 
grated Pedersen conductivity was 23.8 mho. On the down leg, 

- Q 

the Hall conductivity at the single point measured was 5. 3 x 10 

— O 

mho/meter and the Pedersen conductivity was 5.0 x 10 
mho/meter. 
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Model conductivity profiles for the SNT 7/2 flight were 
calculated using exactly the same technique that was used 
for the 7/1 flight. These profiles are also plotted in Figure 17. 
There are substantial discrepancies between the data and the 
model. These differences are probably due to changes in the 
effective conductivity in the earth-fixed frame caused by the 
presence of neutral winds. 

6. Neutral Winds 

The neutral wind velocity has been deduced from the SNT 7/2 
data using the same method that was used for 7/1. The wind 
velocity is plotted in Figure 18 as a function of altitude. 

The scatter in the results from the various analyses is 
much larger in Figure 18 than in Figure u. This implies that 
the systematic errors are relatively larger in the 7/2 neutral 
wind measurements than in the 7/1 measurements. In particular, 
the data below 95 km shows a much larger wind velocity than 
has previously been reported for this altitude range, which 
could be a manifestation of the "settling down" effects mentioned 
above. 

These measurements indicate that the winds were blowing 
in a north to northwest direction on the upleg, with northward 
velocities of about 300-500 meters/ second at all altitudes. 

The westward velocity is small or zero below 115 km and increases 
to around 350 meters/ second above that. The downleg wind was 
southwestward, at about 400 meters/second. Since the trajec- 
tory was northward, the wind configuration was quite remarkable. 
Its implications are discussed in the next section. 
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E. Discussion 

It is profitless and impossible to discuss the electric field, 
current density, and conductivity results from the flight of SNT 
7/1 further until the data from the other experiments has been 
analyzed. 

A generally similar statement can be made about the results 
from the SNT 7/2 flight. However, some comments can be 
made about the implications of the current density measurements. 
As reported, the rocket measurements did not agree with the 
ground magnetometer results. The measurement of strong 
northward and southward components on the up and downlegs 
respectively, when contrasted with the ground measurement of 
no east-west magnetic disturbances, suggest that the rocket 
passed through a current configuration similar to that proposed 
by BostrBm (1968) and seen previously by Park and Cloutier 
(1971), Cloutier et al . (1973) and possibly also on the Javelin 
flight discussed above. This configuration is one where the 
north- south currents flow only across the regions of enhanced 
conductivity associated with arcs, and are connected to the deep 
magnetosphere by east-west oriented Birkeland sheet currents 
of substantial extent, thus forming solenoid-like configurations 
which cause little net external field change. 

The measurement of eastward current on the upleg and net 
westward current from the ground implies that there were 
at least two electrojets present--one eastward and a stronger 
one westward. The fact that the aurora displayed several 
parallel arcs as shown in Figure 3 suggests the hypothesis that 
the total current configuration was composed of several alter- 
nating east and west electrojects associated with several 
alternating up and down Birkeland currents. A similar con- 
figuration was reported by Cloutier etal. (1973). A more 
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detailed investigation of this hypothesis will be possible after 
the completion of the analysis of data from the other experiments 
on the rocket. 

It is difficult to discuss isolated neutral wind measurements 
without more extensive information about the surrounding meteor- 
ological circumstances than is available. Nonetheless, there 
are a few statements which can be made. First, the measure- 
ments of neutral wind reported here are substantially higher 
than any measured in the midnight sector by Brekke et al. 
using the Chatanika radar. Since the analyses performed here 
were done in such a way as to insure maximum comparibility, 
this is curious. It should be noted that the radar results 
represent both a temporal and a spatial average of the data. 

Since other point measurements, such as those reported ty 
Smith, Rees, Stoffregen et al . , and Meriwether et al . , also 
give higher velocity results than Brekke et al . , something 
may be wrong with their procedure. In view of the substantial 
altitude variation shown in Figures 28 and 35, the problem may 
be associated with the height averaging required by their method. 

The results reported here also have some implications 
about the forces which drive high latitude, high altitude neutral 
winds. Historically, the main forces that have been included 
in theoretical studies of these winds were solar heating, tides, 
and inertia (Rishbeth, 1972). Recently, wind data from the auroral 
zone and polar cap and from periods of magnetic storms 
(see the papers cited above) has become available, and has 
not agi eed with such models. Two hypotheses have been advanced 
to explain the disagreement. One hypothesis is that the ano- 
malies are due to pressure forces created by heating of the 
ionosphere by the aurora (Stoffregen, 1972). The other hypo- 
thesis is that the discrepancies are due to ion drag forces and 
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electric field (Fedder and Banks, 1972; Meriwether et al.,1973). 
The data shown in Figure 11 and IkB are not in agreement with 
model predictions given by Rishbeth. Presumably, this was 
caused by one or both of these forces. 

The winds measured on the 7/1 flight could have been driven 
by either force, although the large westward component argues 
somewhat in favor of the ion drag model. The winds detected 
on the 7/2 flight could not have been caused by pressure 
gradient forces. The winds shown in Figure 18 are directed 
at the arc, not away from it as predicted by the heating model. 

In fact, the winds are oriented so as to contribute to building 
up the pressure ridge under the arc. This interesting configur- 
ation appears to have come about because the aurora was rela- 
tively steady for 1 i to 2 hours before the flight. This was 
sufficient time, according to the calculations of Fedder and 
Banks, for the highly structured local electric fields to sub- 
stantially alter neutral wind patterns down to as low as 100- 
110 km. 



Experiment Design Details -- Swedish Nike Tomahawk 7 
Instrument 
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IV. 


A. The Probe 

The decision to use probes that were as small as possible, 
which was made after the first flight, represented a substantial 
mechanical design challenge. The problem of designing a 
probe that was small compared to the Debye length, strong 
enough to withstand the stresses of launch, spin-up and 
deployment, and which could be fabricated with the required 
area accuracy was extremely formidable. Compounding the 
problem was the requirement that the detecting surfaces be 
covered with aquadag with the same area accuracy and without 
the aquadag ever being touched. 

The solution ultimately employed on the Swedish rockets 
is shown in Figure 19. The charge collector plates were 
precision ground tungsten plates. These were glued as shown 
to a less precise substrate, and then were coated with aquadag 
by hand painting under a microscope. This approach produced 
the required strength and area accuracy. However, a certain 
amount of tungsten was left bare of aquadag. This was the 
major compromise in the design. 

B. The Electronics 

The major electronic design challenge posed try this 
instrument was the necessity of measuring a difference between 

4 

two probe currents that might be as small as one part in 10 . 

A block diagram of the resulting circuit is shown in Figure 20. 
The input stage in the probe circuit was, of course, the most 
critical. The input operational amplifiers had to meet quite 



severe specifications. The common mode rejection ratio 

had to be greater than 80 db at 4 Hz, and the bias current 

had to be less than the smallest difference current that it was 

-13 

desired to detect, 5 x 10 amperes. In addition, the offset 

voltage had to be trimmable without increasing the bias current 

significantly. Finally, the slew rate had to be fast enough 

and the parasitic capacitance small enough so that the output 

4 

settled to within one part in 10 of its final value faster than 
the telemetry system could respond using the maximum feed- 
back resistance. This was not easy to do. The amplifiers 
finally used. Intech A-125, met the specifications, but had 
the disadvantage of being quite large (see Figure 22) and 
expensive. 

Obtaining the required dynamic range at the output of 

these amplifiers was also a problem. The possibility of 

using a logarithmic circuit with a diode feedback element was 

explored, but it proved impossible to obtain matched elements 

which were matched with sufficient accuracy. The solution 

used was a set of 8 feedback resistors, which were switched 

in and out of use by mechanical relays driven by an automatic 

gain control circuit. The most important specification on 

the resistors (RO-7 A & B in Figure 21) was that they be 

4 

matched to within one part in 10 . Matched resistor pairs 
of this accuracy are available commercially up to 20 MU. 

Use of a high precision 5 place digital ohmmeter made it 
possible to trim the resistor pairs less than 10 MQ to within 

5 

5 parts in 10 . 

The requirements of the circuit also placed special 
constraints on the relays used to switch in the resistors. 
Parasitic capacitances had to be as small as possible, leakage 
resistance had to be greater than 3 x 10 11 Q and the close 
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and bounce time had to be less than 3 milliseconds. As it 
turned out, only one available relay with the necessary speed 
and lack of bounce also had the high leakage resistance required. 

Also on the input cards were comparators which took the 
difference between input currents and which subtracted the 
sweep voltage. 

In addition to the split probes, the instrument had two 
other probes attached to it. First, there was the so-called 
An experiment shown in Figure 21. This device is the work 
of another experimentor in the Space Sciences Laboratory 
and is outside the scope of this thesis. Second, there was the 

probe. This was an unloaded axis -symmetric probe used 
to provide a stable reference voltage for the sweep. As can 
be seen in the interconnection diagram (Figure 20), this voltage 
was sensed, low pass filtered, and added to the sweep voltage 
on the Gain State Card (Figure 23). 

Also on the Gain State Card was the automatic gain control 
circuitry which sensed the output voltage of the single plate 
current comparator and adjusted the feedback resistance up 
or down if the voltage was outside set limits. Note that 
different gain comparators were sensed depending on whether 
the current was positive or negative. 

The logic of the gain control circuits was arranged so 
that the gain could change one step per telemetry word only. 

It could change every word as long as it continued to increase 
or decrease monotonically . To change direction, it had to 
wait two words unless it was in the highest or lowest states 
possible. 

The last circuit on the Gain State Card was a sample and 
hold and a multiplexer which sequenced the data into the analog 
to digital converter shown in Figure 24> The multiplexer 
and sample and hold are driven by the main programmer card. 
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shown in Figure 25, which also drives the A-D converter, the 
pulse code modulation encoder and readout, the sweep generator, 
and the parallel probe spin motor controller. 

The readout card, shown in Figure 26, takes the various 
quantities passed to it by the ADC and programmer cards and 
assembles them into a 16 kHz pulse code modulated (PCM) 
signal for transmission to the telemetiy system of the rocket. 
The PCM signal was transmitted to ground on FM-FM telemetry 
subcarrier channel H. The main reason for this complex 
procedure was to achieve large dynamic range in the difference 
current channels without absorbing too marry telemetry channels. 
This was achieved by using a 12 bit PCM format (see Figure 28 ) 
and by changing gain every two frames if necessary. The 6 
millisecond frame interval of the PCM telemetry, which was 
set by the bandwidth of Channel H, in turn established the 
basic performance requirements imposed on the input relays 
and operational amplifiers. 

Figure 27 shows the details of the last card in the box, 
the motor controller card. This card accomplished two major 
functions. First, a suitably programmed digital to analog 
converter generated the sweep. Second, drive circuitry 
controlled a stepping motor which rapidly reversed the orien- 
tation of the parallel probe at the start of every fourth sweep. 

This was done to determine the values of any DC offsets which 
the probe had. 

The circuits of the two main pieces of ground support 
equipment are shown in Figures 29 and 30. 
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TABLE I 


FLIGHT 


PROBE GEOMETRY 


PROBE 

CONFIGURATION 


PROBE DIMENSIONS BOOM REMARKS 

LENGTH 


Javelin 8. 56 

Cylinder 

Diameter = 3 mm 
Length = 9. 5 mm 
& two 14 mm long 
guard cylinders 

m 


Swedish Nike 
Tomahawk 7/1 & 2 

Rectangular Prism 

Width = 5 mm 
Length = 10 cm 
Thickness =3.8 mm 
& two 2. 54 cm guard 
sections 

m 

Parallel Probe 
motor mounted 


CO 



TABLE II 


SWEEP CHARACTERISTICS 


FLIGHT PROBE WAVEFORM PERIOD (S) VOLTAGES 


Javelin 8. 56 


Swedish Hike 
Tomahawk 7/1 


Swedish Nike 
Tomahawk 7/2 


Parallel 


1 second 


12 volts 


Perpendicular 

Parallel 

Perpendicular 


Parallel 

Perpendicular 


t r 


0. 1 seconds 


12 volts 


Stepped Ramp 

rt 


6. 3 seconds overall 
3 step pattern 
step widths = . 79 sec. , 
4. 7 sec. , .79 sec . , 


2 volt ramp 
centered on + 1. 7 
steps are 3 volts 
(sweep goes from 
-2.3 to +5. 7 volts) 


6.3 sec. overall 
, 2 step pattern 
step widths = 

' . 79 sec. , 5. 5 sec. 


2 volt ramp 
step is 2. 5 volts 
(sweep goes from 
-2.3 to +2. 0 volts) 


CO 

cn 



Figure Captions 


Figure 1 : 


Figure 2: 


Figure 3: 


Figure 4: 


Figure 5: 


Figure 6: 


Artist's conception of the SNT 7 payload just 
after deployment. (Figure courtesy of SAAB). 

Magnetogram from ESRANGE, Kiruna, Sweden 
for the night of 7-8 February, 1972. Dashed line 
shows time of launch of SNT 7/1. 


M All-Sky M Camera photograph taken from 
JukkasjKrvi, Sweden at 00 56 UT on 9 February 
1973, showing auroral arcs penetrated by rocket 
SNT 7/2. 

IVkgnetogram from ESRANGE, Kiruna, Sweden 
for the night of 8-9 February, 1973. Dsshed 
line shows time of launch of SNT 7/2. Figure 1 
was taken 2 minutes later. 

Plasma density measured by the Uppsala 
Ionospheric Observatory experiments on SNT 7/l. 
(Data supplied by B. Holback, UJO) 


Perpendicular components of the electric field 
measured by the split Langmuir probe on SNT 7/1 
deduced using various assumptions on the parallel 
ion flow velocity in the earth -fixed frame of 
reference. 



Figure 7: 
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Perpendicular electric field components in the 
earth fixed geomagnetic north and west directions, 
and parallel ion flow velocity in the earth-fixed 
frame. Four sets of curves are plotted for the 
electric field, showing field values deduced for 
three assumptions of parallel flow and field 
values deduced from ion flow data "corrected” 
by subtracting Vm^/m^ times the apparent 
electron flow data. (SNT 7/1 data) 

Figure 8: Perpendicular electric field components in the 

earth-fixed geomagnetic north and west directions, 
measured by the split Langmuir probe and by the 
double probe detectors on SNT 7/1; and parallel 
ion flow velocity in the earth-fixed frame 
measured by the split Langmuir probe. The split 
Langmuir probe analysis was done assuming the 
values of v^ as shown. The double probe analysis 
was done assuming the parallel component of 
the electric field is zero. 

Figure 9: Perpendicular current density components 

in the earth-fixed geomagnetic north and west 
directions deduced from the ion flow measured 
by the split Langmuir probe and from the electron 
flow deduced from the electric field measured 
by the double probe detectors on SNT 7/1. 

The consequences of making alternative assumptions 
about parallel ion flow and wake corrections are 
shown for selected altitudes. 

Figure 10: Hall and Pedersen conductivity deduced from the 

current density shown in Figure 26 and the electric 
field measured by the double probe detectors on 
SNT 7/1. The solid lines are the conductivity 
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Figure 11: 

Figure 12; 
Figure 13: 
Figure 14: 

Figure 15: 
Figure 16: 
Figure !7: 
Figure 18: 


profiles predicted by the theoretical model of 
Brekke et al . , (1973) using the density data shown 
in Figure 22. Alternative assumption results 
are shown as above. 

Perpendicular neutral wind components in the earth- 
fixed geomagnetic north and west directions 
deduced from the ion flow velocity measured by 
by the split Langmuir prob^ and the electric 
field measured by the double probe detectors on 
SNT 7/1. Alternative assumption results are 
shown as above. 

SNT 7/2 plasma density. Details are the same 
as Figure 22. (Data supplied by B, Holbaek, UJO) 
SNT 7/2 parallel flow precession analysis. 

Details are the same as Figure 23. 

SNT 7/2 perpendicular electric field components, 
for five different analyses. Details are the'same 
as Figure 24 for four of the sets of curves 
shown. In addition, a set of curves is plotted 
showing the field values deduced using the parallel 
detector data. 

SNT 7/2 perpendicular electric field components. 
Details are the same as Figure 25. 

SNT 7/2 perpendicular current density components. 
Details are the same as Figure 26. 

SNT 7/2 Hall and Pedersen conductivities. Details 
are the same as Figure 27. 

SNT 7/2 perpendicular neutral wind components. 
Details are the same as Figure 28. 
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Figure 19"; 
Figure 20: 
Figure 21 : 
Figure 22: 
Figure 23: 
Figure 24: 
Figure 25: 
Figure 26: 
Figure 27: 

Figure 28: 
Figure 29: 
Figure 30: 


The SNT 7 probe design. 

SNT 7 experiment: interconnection diagram. 

SNT 7 experiment: An circuit design 

SNT 7 experiment: input card circuit diagram. 

SNT 7 experiment: Gain state card circuit diagram. 

SNT 7 experiment: A-D converter card circuit diagram. 
SNT 7 experiment: Programmer card circuit diagram. 
SNT 7 experiment: readout card circuit diagram. 

SNT 7 experiment: motor controller card circuit 
diagram. 

SNT 7 experiment: telemetry format. 

SNT 7 experiment: Check out box circuit diagram. 

SNT 7 experiment: Calibration box circuit diagram. 
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